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Abstract 

We study quantum filters that are driven by basic quantum noises 
and construct classical versions. Our approach is based on exploiting 
the quantum markovian component of the observation and measure- 
ment processes of the filters. This approach leads in a natural way the 
classical versions for a class of quantum filters. We consider quantum 
white noises derived from Wiener and Poisson processes that drive the 
signal and measurement processes and derive the recursive filtering 
equations using classical machinery. 

1 Introduction 

Interest in quantum stochastic filters is fueled by experiments in quantum 
feedback control [3] techniques. There are also theoretical treatments of 
such processes and one such approach is based on the framework of Hudson- 
Parthasarathy version of quantum stochastic calculus [10] . For an excellent 
introduction to quantum filters and the theoretical modeling based on quan- 
tum stochastic calculus the reader is referred to the refernce [8j . The authors 
have derived the recursive relations for the quantum stochastic filters driven 
by quantum white noises for the case of homodyne detection and photon 
counting experiments. Even though the underlying noises are quantum in 
nature the resulting filters are classical stochastic processes. Intrigued by 
this result we set to explore the possibility of deriving the filtering equations 
within the framework of classical stochastic methods. Our hope is that the 
rich literature available in classical stochastic processes can be exploited to 
gather further insights into quatum stochastic filters. As our aim is to es- 
tablish the applicability of the classical techniques we will focus on deriving 
the equations presented in [S] and work with the physical model and the no- 
tations used there in. For the classical techniques we will adapt the methods 
presented in [T], which is a freely downloadable book on the introduction to 
stochastic calculus and filetering, to derive the equations for homodyne de- 
tection. For the case of photon counting problem we will use the techniques 
from the reference jllj . 



Our work is motivated by the research that relates classical and quantum 
stochastic processes and the influence each field has on the other. The rela- 
tion between non-commutative and classical probabilities has been explored 
since 1970 O 0] and continues to grow at a steady pace [E]. Quantum 
probability which is a non-commutative probability defined on a von Neu- 
mann alegbra [9]. It is well known that we can construct classical versions 
from a quantum probability space using spectral theorem if the underlying 
*-algebra is commutative [8j. There are ways to realize the classical versions 
even if the underlying algebra is a non-commutative von Neumann alge- 
bra. As has been established in [3] a restricted class of quantum markovian 
processes can be realized as classical versions. Our report is organized as 
follows: first we describe the white noise processes that make the physical 
model in detail and establish them as quantum markovian in the sense of 
[3]. Then we build the corresponding classical versions and prove that they 
have the same covariance as their quantum counterparts by adapting the 
techniques from [5] . Finally we derive Martingales on a classical probability 
space for the measurement processes and express them as integral of some 
Brownian motions. This will establish the classical counter part for our 
quantum filtering problem. [1] 

2 Quantum stochastic calculus and filtering 

We start with the physical model described in [8] for the quantum filtering 
problem that is applicable to quantum optics. An example of such a system 
would be an atom placed in an electromagnetic field. The model consists 
of an intial system defined on a finite dimensional Hilbert space that 
interacts with an environment described by a symmetric boson fock space 
r = r5(L^([0, T)). Let us work with the quantum probablity space (.4, P) 
where ^ is a von Neumann algebra defined as (8) where W = i3(r), 

the von Neumann algebra generated by the bounded operators on F. We let 
P as a normal state defined as (p. In this work we will assume the state 
(j) to be the fock vacuum state and all the operators to be bounded. Thus 
the state of our von Neumann algebra is faithful normal. 

Accordingly, the quantum stochastic differential equations describing the 
filetering problem are 



djt{X) = jt{CLAX))dt + 3t{[L\X])dAt+ jt{[X,L])dAu (2.1) 



where X is an obeservable that belongs to the bounded set of operators 
of J^Q and jt{X) = Ui{X®)Ut describes the evolution of the observable in 



A 



dY^ = jt{L + L*)dt + dAt + dAl; 
dAt + jt{L)dA; + jt{L*)dAt + jt{L*L)dt; 



(2.2) 
(2.3) 
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Heisenberg picture. The unitary operator Ut satisfies the quantum stochastic 
differential equation 

dUt = {LdA; - L*dAt - ^L*Ldt - iHdtjdUt, Uq = I (2.4) 

In the above, equation ()2.ip represents the signal process, ()2.2p repre- 
sents the homodyne detection, and (j2.3p represents the photon counting 
measurement prcoesses. The filetering problem is stated as the derivation 
of conditional law TTt{X) = ¥(jt{X)\'3^t) where is the von Neumann alge- 
bra generated by the observation process which is commutative and thus a 
classical process. 

Let us consider the processes, Zt{f) = atAt{f) + PtA*{f) and Wt{f) = 
Zt{f)+^t^t{f) where At, A^, and Af are the creation, annihilation, and con- 
servation or gauge processes of the fock space, at, fit, and 7^ are time depen- 
dent complex variables, and / € T . Let us call them as ^/K-Gaussian (Non- 
commutative) and c/K-Martingale processes. These are non-commutative 
processes due to the arbitrary time dependent coefficients and we will see 
the conditions under which they posses classical versions. These processes 
satisfy the differential equations 

dZt = atdAt + PtdA*t (2.5) 

dWt = dZt + itdht (2.6) 

Let us take a closer look at the structure of the Fock space operators. In 
the context of quantum optics they are defined as and A*^^^ where g is 
the charecterestic function 1[q^^] satisfying CCR on the exponential domain, 
and on the finite particle vectors with covariance [10] 

[Ag{t),Al{s)]=min{s,t) (2.7) 

The gauge process which is essentially the Poisson process is defined as 

(AJ)(r) = Aft(r)/(t) = |rn[0,t]|/(r), / E T, r G f), t G [0, T] (2.8) 

It is clear that equations (j2.5p an(i (|2.6p form the components of the signal and 
measurement processes. Following [4J we can now define the ^/K-analogues 
of classical Gaussian processes based on their covariance. 

Definition 1. The ^-Gaussian process X't^^'^ with covariance min(s,t) is 
called a ^-Brownian motion. 



As we mentioned earlier we adopt the notion of quantum markovian, 
quantum martingale, and independence defined in jlj. 
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Definition 2. Let (^,P) be a quantum probability space and {X)t^T is a 
stochastic process on P). Let 

£/^^-=vN{Xu\u<t) (Z£^ (2.9) 

^j^] := vN{Xt) C £^ (2.10) 

The process Xt^x is called quantum markovian if the following holds: Vs, t G 
T with s < i 

P[X|A]] C^[,], VXG^[i] (2.11) 

The conditional expectation is a map from A to Agy In our context 

this becomes, at Hilbert space level, a projection from Jif® T to J^® F^j 
(see [1] for details) . 

Definition 3. Let (i2/,P) be a quantum probability space, a process [X)t^T 
defined on it is a called a martingale if 

P[Xt|=<]] =X„ Vs<t (2.12) 

The following proposition will be helpful in establishing the martingale 
property of ^-Gaussian processes. 

Proposition 1. (Ref f^) A ^ -Gaussian process is a martingale if and 
only if Pfs] = /sVs < t which is the case iff c(s,t) = c(s,s) Vs < t. 

Proposition 2. (Ref I4j) ^ -Brownian motion and jV -Martingale pro- 
cesses are quantum markovian as well as quantum martingale. 

Proof. A ^-Gaussian process with covariance c is quantum markovian iff 
V s,u,t G r, such that s < n < t, 

c(t, s)c{u, u) = c{t, u)c{u, s) (2-13) 

Thus cy4^-Brownian motion processes is quantum markovian. See the proof 
of theorem 4.9 in [Hj for the validity of equation (10). 
We have c(t,s) = min(s,t) = s = c(s,s), the martingale property is estab- 
lished. 

The gauge process being the Poisson process the markovian and martingale 
properties of ^-Martingale follow immediately. □ 

Having established the quantum markovianity of the ^-Gaussian and 
^-Martingale processes we will build the classical version that preserves 
the covariance. 
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Definition 4. Let {Xt)t^T be a stochastic process on a quantum probability 
space {£/,¥). We call a process (Vj)tgT taking values in M and defined on 
a classical probability space {Q,'^ , ^) a classical version of {Xt)t£T 

if all 

time ordered moments coincide. That is Vn E Wti, . . . ,tn S T with 
ti < ■ ■ ■ < tn, and all bounded Borel functions hi, hn on M the following 
equality holds: 

F[hi{Xt,). . . KiXtJ] = [ hi{Vt,{u^)) . . . K{VtM)d^{Lo) (2.14) 

Jn 

Remark 1. Our interest in time ordered moments is due to its applicability 
in covariance which is needed to characterize Gaussian processes. From the 
results on time ordered moments we infer covraince on two different epochs 
are preserved irrespective of their order. 

Our classical probability space is the real line ^ and the sigma field is 
the sigma field generated by the borel sets. We define the probabilty measure 
and the classical version Vt as ^{Vt £ B) = P(lB(Xt)), where B e . Let 
us define the transition operators for the quantum markov processes that is 
the same in the classical context. Then we realize the classical version using 
the transition operator, a la Kolmogorov construction. This a widely used 
mathematical technique as can be seen from [1], [2], and |12j . 

Definition 5. Let {Xt)t£T as defined above and is self-adjoint \/t € T. Let 
spect{Xt) and vt are the spectrum and distribution of Xt respectively. Let 
us define the space 

L^{Xt) = vN{Xt) = L°°{spect{Xt), vt). (2.15) 

and the transition operator 

J^,t : L°°{Xt) — > L°°{X,), with s<t. (2.16) 

defined as 

p[/i(xoK.]] = p[/i(Xt)K[,]] = {jes,th){Xs). (2.17) 

Theorem 1. The -Gaussian process X"^*^ has a classical version Vt- In 
particular the covraince is preserved for the Vt Gaussian process and hence 
is a classical Brownian mtion. 

Proof. Following the line of reasoning in [3] theorem 4.4 we get 

F[hi{Xt,), hniXtJ] = n{hi.^t,M{h2-^t,M{hz. . . . ))){Xt,)\ (2.18) 
It follows that for arbitrary ti and t2 we have 
¥[hi{Xt,),h2{Xt2)] = ^[hi{Vt,),h2{Vt,)]=c{ti,t2) = min{tiM) (2.19) 



In essense we have replaced X" that has non-commutative components 
with a single random variable such that the first order moments coincide. 
The resulting stochastic process is commutative on a classical probability 
space. □ 
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3 The filtering problem 



The objective is to find a least mean square estmate of a system of observable 
X £ ^ at time t, given the observations Yt up to this time, that is to obtain 
a recursive relation 

7rt{X)=F{jt{XWt) (3.1) 

where S^t = vN{Ys : < s < t} is the commutative von Neumann algebra 
generated by the observation process. Let us now adopt the definition of 
conditional expectation from 

Definition 6. (Conditional expectation). Let be a quantum prob- 

ability space and let ^ C he a commutative von Neumann subalgebra. 
The conditional expectation P(.[^) is a map from ,yi^' onto ^yV such that 
the following holds: 

F{.\.yy) -.^'^.yy, 

F{F{B\^)A) = F{BA), MAe J^.B e JV' (3.2) 
Here jV' is the commutant of . 

Remark 2. ^ is the commutative von Neumann algebra generated by the 
observation process. That is, an algebra of compatible observables. The 
domain of the conditional map has to be the commutant of in order to 
make sure that compatible observables are present in the same realization. 
If in addition B = B* then B is also in the commutative algebra and hence in 
the classical probability space by spectral theorem. The self-nondemolition 
property is established easily as in [8]. 

We can now define all the objects of ()2.ip . (12. 2p . arad ()2.3|) to the corre- 
sponding classical objects and derive the filtering equations using the classi- 
cal results as in [Ij and for the photon counting problem based on [11]. Let 
us rewrite them using ()2.6p ared (j2.6p as 

djt{X) = jt{CL,H{X))dt + dZt; (3.3) 

dVt^ = jtiL + L*)dt + dZf, (3.4) 

dY/" =jt{L*L)dt + dWt; (3.5) 

Theorem 2. (Kushner-Stratonowich, ref II]) Let X and Y are the signal and 
observation processes satisfying the following classical stochastic differential 
equations 

dXt = fit, Xt)dt + a{t, X)dVt (3.6) 

dYt = h{t, Xt)dt + dWt, Yq = (3.7) 

where Vt and Wt are two independent Brownian motions. Then the filtering 
equation is given by 

d7Tt{<P) = 7:t{As{<P)) + (vrt(/i^,/.) - 7rt{h^)7Ttimdyt " Mh)dt) (3.8) 
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Theorem 3. (ref f77] /} Let Xt,t > is a semimartingale and Yt,t > is a 
counting process such that 

Xt = Xo+ I Hsds + Zt (3.9) 
Jo 

where Zt is a ^t-'^O'fii'ngale, Ht is a hounded ^t-p^ogressive process, E[sups<t\XsW < 
oo for every t > and Xq is a square integrable randaom variable. 

Yt = YQ+ [ hsds + Mt (3.10) 
Jo 

where Mt is a .^t-'m'^ftingale with mean and a{Mu — Mt; u > t) is indepen- 
dent of a{Yu, hu] u >t) and ht = h{Xt) is a positive bounded ^t-progressive 
process such that E h'^ds < oo,\/t. Then the filtering equation is given by 

7T{Xt)=7T{Xo)+ l\{Xs)ds+ [\7T{hs)]-'[7T{X,--hs)-7T{X,^)7T{hs)]dms, 

Jo Jo 

(3.11) 

where ^ 

mt = Yt [ 7r{hs)ds (3.12) 
Jo 

Now we are ready to derive the main result of classical versions for 
quantum filters. 

Theorem 4. A classical version exists for the quantum filtering problem 
defined by the equations (j3.3p . (|3.4|) . ared (j3.5|) . 

Proof. Let us apply spectral theorem on the commutative von Neumann 
algebra (^,P) of the observation process (Ref [1^ Proposition 1.18.1). The 
abelian ^ is *-isomorphic to a W* algebra L°°(J7, J^, //). Let us outline some 
of the main steps of the proof. 

The linear functional P defines a Radon measure such that 

P(a) = / a{t)dfi Va G ^ (3.13) 
Jn 

The mapping $ : a — > a{t) is a *-homomorphism such that 

¥{ab)= [ a{t)b{t)dii ^a,be'¥ (3.14) 
Jo. 

Therefore the target random variables a(t) are square integrable. 

We then apply Radon-Nikodyn theorem to get a probability measure (p 
that is absolutely continuous with respect to /i. In the vacuum state the 
quantum Martingale component of the observation processes are mapped 
to classical square integrable Martingales, a continuous Martingale in the 
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case of Homodyne detection and a cadlag Martingale in the case of photon 
counting, with respect to the natural filtration, that can be represented as 
integral with respect to Brownian motions. 

On the same probability space using Thoerem 1 we can replace the non- 
commutative process Zt with its classical version Zt which is a Brownian 
motion. 

So the classical version of the signal and observation processes become, 

d3t{X) = Y.^{i{^l,h{X))) + dZt (3.15) 

dW^ = jt{L + L*)dt + dVt (3.16) 

= jt{L*L)dt + dilt (3.17) 

where Z, V , U are all independent Brownian motions. We get by applying 
Theorem 2 for the Homodyne detection and Theorem 3 for photon counting 
respectively, 

dntiX) = ^t(^L,j^)(it+(^t(L*l+XL)-7rt(L*+L)7rt(X))(d^-^t(L*+L)dt) 

(3.18) 

dTTtiX) = M^L,H)dt + C'^^lfJ:^ - MX)){d&t - ML*L)dt) (3.19) 

□ 
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